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In a recent work on low energy pion-nucleon scattering, instead of using chiral perturbation theory 
(ChPT) amplitude, we started from a pion-nucleon soft-pion result and used elastic unitarity directly 
as a dynamical constraint to construct first-order unitarity corrected amplitudes. The resulting am- 
plitudes are crossing symmetric but, as the ChPT ones, satisfy only approximate unitarity relation. 
In the present work, we reconsider our approach and we apply the inverse amplitude method (IAM) 
in order to access the energy resonance region. We present the resulting S- and P-wave phase shifts 
£h , that are shown to be in qualitative agreement with experimental data. 

a , 

O ! I. INTRODUCTION 

cn . 

■ Even though Quantum Chromodynamics (QCD) is considered nowadays the theory of strong interactions, its 
^ ' application to low energy hadron physics is still far from a solved problem in physics. A great theoretical improvement 
was made by means of the method of Chiral Perturbation Theory (ChPT) [1], which is an effective theory derived 
from the basis of QCD. The method consists of writing down chiral Lagrangians for the physical processes and uses 
the conventional technique of the field theory for the calculations. It is a quite successful method when applied to 
^— *j • meson processes. 

| In order to deal with baryons in the ChPT approach, the theory faced problems related to fixing the scale for 
momenta and quark mass expansion. This led eventually to a method known as Heavy Baryon Chiral Perturba- 
tion Theory (HBChPT) [2], which has been applied to describe pion-nucleon scattering from a complete effective 
Lagrangian calculated up to third order in small momenta [3]. More recently the explicit degrees of freedom corre- 
sponding to the A and N* resonances have also been considered within HBChPT [4]. One peculiar feature of ChPT 
Oh, is that it leads to partial waves satisfying only approximate elastic unitarity relation, since, for instance, the leading 
amplitude is a real function of energy in the physical region. 

Unitarization methods have been applied to pion-nucleon scattering in the literature for a long time [5]. An 
interesting approach to implement unitarity in HBChPT amplitude is the Inverse Amplitude Method (IAM) [6]. This 
is a sort of N/D method that takes the leading contribution to the partial wave amplitude as the numerator N and 
includes the imaginary part that comes from loop corrections in the denominator D. The expansion of D yields the 
same structure of the original amplitude up to the order considered, plus corrections of higher orders, but is fully 
unitary. IAM partial wave amplitudes for pion-nucleon scattering have been constructed and fitted to the experimental 
phase shifts by fixing nine free parameters [7] . 

On the other hand, instead of working in an effective theory framework, one can use the hard-meson method of 
current algebra. In fact, current algebra approach was abandoned by several reasons but, as it is based on chiral 
symmetry, one expects that it gives equivalent results for meson processes. In fact, we have compared these two 
methods and concluded that quasi-unitarized current algebra amplitude is equivalent to ChPT one-loop calculation 
in the case of kaon-pion scattering and equivalent to one- and two-loop ChPT calculations for pion-pion scattering 
[8]. In particular, both methods lead to amplitudes satisfying the same approximate unitarity relations. 

Two years ago, we developed an alternative procedure to unitarize a current algebra pion-nucleon scattering am- 
plitude [9]. The unitarity correction was built as follows. We started from a soft-pion amplitude [10] reproducing 
Weinberg prediction for S-wave scattering lengths [11] and constructed auxiliary functions in order to respect approx- 
imate unitarity. Considering known the imaginary parts of partial wave amplitudes, we used the dispersion relation 
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technique to arrive at a quasi-unitarized amplitude written in terms of known functions and a two parameter polyno- 
mial part. Furthermore, we imposed exact crossing relations to Dirac amplitudes resulting very constrained partial 
waves. By this method, even in the absence of free parameters, the violation of unitarity of the resulting amplitudes 
seems to be very important mainly in the resonance region. 

The motivation of the present work is to go beyond threshold by constructing partial wave amplitudes that respect 
the unitarity relation exactly. In the present paper, we reconsider the unitarity corrected amplitudes described above 
in order to construct IAM modified amplitudes. In the approximation used here, each partial wave requires a maximum 
of two parameters. 

In the next section we present a summary of the basic formalism. In section III we present the procedure of Ref. 
[9] for obtaining quasi- unitary corrected amplitudes. In section IV we present the IAM unitarization procedure and 
the results. 



II. BASIC FORMALISM 

We consider the reaction ir a (k) + N(p) — > 7r b (fc') + N(p'), which is described by the T-matrix amplitude 

T a \p',k'; Pl k) - -A ab (.s,t,u) + + B a \s,t,u), (1) 

with s = (p + k) 2 , t = (k' - k) 2 and u = (p - k') 2 . 

In order to specify the various charge states, the invariant amplitudes A and B are decomposed as 

A ab = S ab A + + l[ r 6 ;T a] A - ; 

B ab = S ab B+ + ^[T b ,T a ]B-. (2) 

They exhibit the following symmetry properties under crossing: 

A ± (s,t,u) = ±A ± (u,t,s), B ± (s,t,u) = =f B ± (u,t,s); (3) 
the amplitudes corresponding to definite isospin 1/2 and 3/2 are given by 

A 1/2 = A+ + 2A-, A 3/2 =A+-A-, 

and similarly for B 1/ / 2 , -B3/2 • 

We work in the center of mass system, so that the four momenta arc defined as 

k = 0,w), k' = {k',w), p= {-k,E), and p' = ( -k',E), 

with \k\ = \k'\, w = Vk 2 + m 2 and E = Vk 2 +M 2 , M = .938 GeV and m = .140 GeV being the nucleon 
and the pion mass, respectively. The total energy and scattering angle are given, respectively, by W = E + w and 
\k\ 2 cos 9 = k-k', thus, in terms of these quantities, one has s = W 2 , t = —2 \k\ 2 (1 — cos8) and u = 2M 2 + 2m 2 — s — t. 
For each isospin I the Pauli amplitudes are 

Fu(a,cos0) - ^^-[ Ai(8,coaO) + (W-M)Bi(8,cos0)], 

F 2I (s,co S 6) = ^^-[-A i (8,cob0) + {W + M)B i {8,cob0)]. (4) 

Partial wave amplitudes are defined as 

ff t {s) = F 1It (s) +F 2ie ±i(s), where 

1 f +1 1 3 

F iid s ) = ^J F u (s,x)P e (x) dx, for i = 1,2 and / = -, -. (5) 

For elastic scattering we have Im/^(s) = \k\ \ff t {s)\ 2 , which may be solved yielding 



fh(s) = ^re iS '^ sinM*). (6) 
\k\ 



where Sie(s) are real phase shifts. 
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III. UNITARIZATION PROCEDURE 



Let us recall the main points in our unitarization procedure. First, we consider an amplitude reproducing S-wave 
scattering lengths predicted by Weinberg. The low energy amplitude we start from is the soft-pion limit obtained 
by Osypowsky [10] using the Ward identity technique. The final expression for the four point function, related to 
pion-nucleon scattering, is written in terms of form factors and propagators. By estimating the contribution of each 
term at threshold one is led to assume that the low energy amplitudes are 

A ca- AV ^ B ca- = L±gv ; A ca+ = B ca+ = 0, (7) 



8M/ 2V n 2f 2 

where / = .094 GeV is the pion decay constant, [i v ~ 3.7 and the superscript ca stands for current algebra. These 
amplitudes lead to the well known Weinberg prediction for S-wave scattering lengths, namely a$ = and 
= 0.077 m _1 , to be compared with the experimental values = — 0.015±0.015 m" 1 and = 0.097±0.003 mT 1 . 
We implemented unitarity for the low-energy pion-nucleon amplitude A ca ~ and B ca ~ . From our previous analysis 
on meson scattering, we conjecture that the corrected amplitudes must satisfy 

A e (s) ~ Af 1 (s) + Af\s) + 0(e 2 ), for s ~ (m n + m N ) 2 , and the same for B e ; 

A^ is a complex function and e is a small parameter characterizing the corrections. 
For each isospin channel, soft-pion amplitudes are obtained from (5) as 

/r + (s) = F ca {s) + F ca (a)) f ca- [s) = pca^ 
f ca- {s) = pca {s) + F ca (s); f ca + {s) = pea ^ 

where f£ a ± follow from A ca and B ca . It is evident that Im ,/ £ ca± (s) = 0, but this situation is inherent to 
soft-meson calculation. 

We construct quasi-unitarized amplitudes by requiring that the corrections F$, for s > (M + m) 2 , have the 
following imaginary parts: 

Im = \k\(Fr Q 2 ( S ) + 2Fft{8)Fft(8)), 

ImF^(s) = \k\(Fi« 2 (s) +2F^l{s)), 

lmF^(s) = |fc|^ c f( S ), Im F${8) = \k\F^ 2 {s). 

It guarantees that the quasi-unitarized partial waves ff — f c t a± + f^ ± satisfy 

Im/± = |fc| |/±«f. (8) 

Expressing functions A and B in terms of Pauli amplitudes, we construct the auxiliary functions 

A(s, cos8 s ) = - [ai {s)S(s) + a 2 (s)V(s) + 3 cos e s a 3 (s)Q(s)] , 

S(s,cos0 s ) = ^[b 1 ( s )S(s)+b 2 (s)V( s )+3cos9 s b 3 ( s )Q( s )], (9) 



where 



Im S(s) = 3*1 A?{s), ImP( S ) = ^C Q(s) = ^A?(s), 





= (W + 


M) ( F™ + 2F™) + (W - M) (F 2 c g + 2Fft) , 
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= (W 2 - 
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In order to avoid kinematical singularities, we write subtracted dispersion relations for S, T> and Q by introducing 
free parameters. At this point the present procedure must be handled in a different manner than done in Ref. [9]. In 
the construction of the quasi-unitarized amplitudes we introduced only two parameters. On the other hand, in the 
present work, one must use up to two parameters for each partial wave, chosen from the set of parameters (Ai, A 2 , A3) 
in the expressions 

S(s) = s 2 A! + A c a (s) G(s), V{s) = s 2 A 2 + B c a G(s), Q(s) = s 2 X 3 + A™{s) G(s), 

with 



n J(. 



x — (M + m) 2 ] [x — (M — m) 2 



(M+m)2 x 4 (x - s) 



dx. 



One may argue that in the definition of the subtracted dispersion relations above, terms in powers of s smaller than 
two are missing. This restricted choice of the number of subtraction constants is based on the fact that in the fitting 
process those extra free parameters do not provide qualitative improvement on the fits. 

As done in our previous work, crossing properties and partial wave total isospin dependence are imposed to the 
corrected amplitudes, by taking 



A^+(s,t,u) = 2A(s,t) + (s^ u), A^~(s,t,u) = A(s,t) - (s 
B (1)+ (.s,t,u) = 2B(s,t) - (s <-> u), B^~(s,t,u) = B{s,t) + (s < 



> u) 
u). 



The corrections to partial wave amplitudes are calculated from (5) using (4) and arc indicated by /^(s). This was 
the final step of the procedure introduced in Ref. [9]. 

IV. INVERSE AMPLITUDE METHOD AND RESULTS 



The quasi-unitarized amplitudes from last section do not obey exact unitarity relation. In order to restore it, we 

fAby 

Wis) 



apply the IAM to the corrected partial waves fa = + f u \ by writing 



///(*) 



1- /#(*)/./??(*) 



In order to fit IAM amplitudes to the experimental results of phase shifts one has to be cautious. As the inelasticity 
for some partial waves becomes important for energies less than 1.5 GeV, one has to introduce some criteria to 
define the validity region of the IAM. In the present paper we decided to fit our results to experimental phase shifts 
corresponding to dimensionless experimental amplitudes f^ p satisfying Im f^ p — \fg^ p \ 2 < 0.1. 

Restricting ourselves to these ranges in energy, we performed the fits of Si 1 , S3 1 , Pi 1 , P3 1 , Pi 3 and P 3 3 partial 
waves from the form f(j above to experimental data, thus fixing the free parameters Ai, A 2 and A 3 independently, 
but only a maximum of two of them for each wave. Their values are presented in the Table. 
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-3.8 


A 2 


24.3 


39.8 




2.3 


-48.0 




A3 




16.3 


-6.9 


-8.5 


-27.0 


-8.3 



TABLE I. Results from fits of IAM modified quasi-unitarized partial waves to experimental data; Ai and A3 are given in 
units of (GeV)~ 5 , while A 2 , in units of (GeV)" 6 . 
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The phase shifts thus obtained are shown in the figures, as functions of the cms energy, in GeV. In the same figures 
we plot the phase shifts obtained with the original quasi- unitarized amplitudes [9]. In that case only Ai and A3 were 
considered, and their values were the same for all partial waves. Concerning our present results, we observe that it is 
possible to fit the phase shifts of our model to the experimental data (Ref. [12]). We also observe that, using the IAM, 
it is not possible to simultaneously fit two waves with the same values of parameters. The qualitative agreement with 
experimental data is an indication that the IAM allows one to access the resonance region of pion-nucleon scattering 
since some low energy chiral amplitude is given. 

We would like to emphasize that in the construction of the quasi-unitarized amplitude there is no commitment on 
the existence of any resonance. They emerge as dynamical consequences of the unitarization procedure. Notice also 
that we do not include the nucleon pole in the low energy amplitude we started from. 

In order to further explore the results that we have obtained, one should construct corrections in the next level of 
unitarity approximation for the partial waves, as outlined in the last equation of Ref. [9]. We believe that the IAM 
applied to this new result will allow one to access the energy resonance region for higher angular momenta partial 
waves. 
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FIG. 1. S-wave phase shifts (in degrees) as functions of cms energy (in GeV); present results (solid) and previous quasi-unitary 
ones (dashed). 
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FIG. 2. Isospin 1/2 P-wave phase shifts (in degrees) as functions of cms energy (in GeV); present results (solid) and previous 
quasi-unitary ones (dashed). 
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FIG. 3. Isospin 3/2 P-wave phase shifts (in degrees) as functions of cms energy (in GeV); present results (solid) and previous 
quasi-unitary ones (dashed). 
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